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FOREWORD 

The  U.  S  Army,  through  the  facilities  of  the  U.  S.  Army  Transpor¬ 
tation  Research  Command,  Fort  Eustis,  Virginia,  has  been  studying  various 
methods  for  extending  the  range  of  Army  aircraft  and,  in  particular,  the 
range  of  the  helicopter.  One  method  being  investigated,  which  indicates 
that  ferry  ranges  of  2,000  miles  or  more  are  attainable,  is  one  in  which 
floating-wing  fuel  tanks  are  attached  to  a  helicopter's  fuselage  through 
a  hinge  connection. 

The  report  presented  in  the  following  pages  is  Volume  3  of  a  five- 
volume  final  report  of  the  analysis  and  wind-tunnel  investigation  of  the 
floating-wing  system.  Volume  3  presents  the  results  of  the  wing  flutter 
analysis  of  the  floating  wing  at  various  fuel  conditions.  The  con¬ 
clusions  contained  herein  are  concurred  in  by  this  Command. 
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SUMMARY 


Anti-symmetric  and  symmetric  flutter  investigations  are  performed  on  the 
Range  Extension  floating  wing  for  various  fuel  loading  conditions  and 
combinations  of  wing  modes .  The  method  employed  in  this  analysis  makes 
use  of  Theodorsen's  oscillating  aerodynamic  loads  for  the  generalized 
excitation.  Mechanical  equations  of  the  system  are  obtained  with  Lagrange’s 
method  using  several  combinations  of  three  wing  modes,  and  two  rigid  fuselage 
modes  as  generalized  coordinates.  Solution  of  the  flutter  determinant  is 
obtained  by  an  IBM  650  computer  program.  Results  of  the  flutter  calcula¬ 
tions  are  presented  as  plots  of  airspeed  against  wing  modes  frequency  ratio, 
from  which  the  flutter  speed  is  determined  by  the  intersection  of  the  actual 
frequency  ratio  with  the  characteristic  curve. 

Results  of  the  analysis  indicate  that  the  flutter  characteristics  of  the 
delta  hinge  wing  arrangement  are  unacceptable.  The  investigation  shows 
that  for  both  the  07.  and  1007.  fuel  configurations  the  flutter  speeds  are 
below  the  cruise  speed  of  80  knots.  It  is  indicated  that  the  low  flutter 
critical  speeds  calculated  here  are  closely  associated  with  the  aircraft 
roll-wing  rigid  body  flap  instability  which  has  been  determined  to  be  an 
aircraft  flight  instability.  The  wing  design  is  being  revised  to  a  straight 
hinge  configuration  in  which  the  necessary  aerodynamic  wing  spring  is  being 
furnished  by  a  differential  flap  system.  This  change,  which  will  eliminate 
the  unusual  flap-pitch  coupling  in  the  oscillatory  airloads,  should  improve 
the  flutter  situation. 

The  effects  of  fuel  slosh  calculated  herein  can  be  evaluated  on  a  rel¬ 
ative  basis.  The  slosh  is  approximately  simulated  by  partial  fuel  con¬ 
ditions  in  which  it  is  assumed  that  the  fuel  is  sloshed  fully  forward  or 
fully  aft  in  the  tanks  so  as  to  cause  maximum  shift  of  the  chordwise  C.G. 
position.  Results  show  that  the  aft  C.G.  position  obtained  from  partial 
fuel  in  all  tanks  lowers  the  basic  flutter  speed  appreciably.  This  can 
be  avoided  by  a  fuel  drainage  schedule  wherein  one  tank  at  a  time  is 
fully  emptied,  tending  to  eliminate  the  critical  aft  C.G.  along  the  entire 
span. 

To  evaluate  the  final  wing  configuration  from  a  flutter  standpoint,  it 
is  recommended  that  a  ground  shake  test  be  performed  to  substantiate  the 
wing  frequency  estimates  used  in  the  analysis.  Then  flight  flutter  tests 
should  be  conducted  starting  with  the  empty  wing  configuration  and  con¬ 
tinuing  with  successively  larger  fuel  quantities  until  the  full  range  of 
fuel  variations  is  explored. 
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I .  INTRODUCTION 


This  flutter  study  is  part  of  an  analytical  and  wind  tunnel  study  being 
conducted  under  the  Reference  1  Transportation  Research  Command  contract. 
The  program  is  aimed  at  the  development  of  a  means  for  helicopter  range 
extension  through  application  of  a  floating  fuel  wing  concept.  Under  an 
earlier  contract,  Reference  2,  Vertol  conducted  an  initial  feasibility 
study  of  the  floating  wing  concept.  Results  of  this  study  are  reported 
in  Reference  3.  The  present  analytical  and  wind  tunnel  investigation  is 
under  the  Reference  1  contract  and  is  based  on  a  Vertol  proposal, 

Reference  4.  Phase  I  of  the  present  contract,  consisting  of  wind  tunnel 
and  mechanical  instability  studies,  was  reported  in  Reference  5  and  6. 

The  present  flutter  analysis  is  a  part  of  Phase  II  of  the  contract. 

The  range  of  present-day  helicopters  with  normal  fuel  load  is  less  than 
400  nautical  miles.  Even  with  additional  internal  tanks,  the  helicopter 
range  is  less  than  1100  miles.  With  floating  wings,  the  range  can  be  ex¬ 
tended  to  as  much  as  2400  miles,  corresponding  to  the  longest  over-water 
distance  on  the  Pacific  Ocean  ferry  route. 

Each  floating  wing  contains  compar tmented  tankage  connected  by  lines  to 
the  helicopter's  main  tank.  The  wing  lift  supports  the  fuel  weight  that 
it  carries,  and  the  helicopter  acts  as  a  tow  to  propel  the  wing  forward. 
Wing  attachment  to  the  helicopter  is  through  a  hinge  so  as  to  eliminate 
the  bending  moments  at  the  fuselage  applied  by  conventional  wings,  thus 
avoiding  the  addition  of  extensive  wing  carry  through  structure  to  the 
helicopter . 

The  hinge  line  is  not  horizontal,  but  is  skewed  aft  as  shown. 


As  fuel  is  consumed  and  the  wings  become  lighter,  they  tend  to  flap  up¬ 
ward  about  the  hinge.  Because  of  the  skewed  orientation,  the  angle  of 
attack  at  any  chord  line  is  reduced  as  the  wings  flap  up,  the  lift  is 
reduced,  and  the  wing  assumes  a  new  mean  position.  Full  span  pilot 
controlled  wing  flaps  are  also  provided,  so  that  the  trim  attitude  of 
the  wing  may  be  adjusted;  these  are  also  used  as  high  lift  devices 
during  the  running  takeoff. 
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The  analysis  performed  herein  is  intended  to  establish  the  flutter  char¬ 
acteristics  of  the  floating  wing  system.  In  normal  flight,  the  surfaces 
of  an  aircraft  are  subject  to  almost  continuous  disturbances  by  random 
aerodynamic  loads.  Gusts  and  pilot  course  corrections,  for  example,  put 
temporary  loads  on  these  surfaces.  As  with  every  elastic  body,  such 
loads  cause  deflections  of  the  aircraft  structure  and,  at  each  application 
of  load,  induce  vibration  in  the  body.  Ordinarily,  these  vibrations  are 
highly  damped  and  rapidly  disappear.  However,  under  certain  conditions 
involving  the  elastic,  inertial  and  aerodynamic  loads,  the  vibrations  are 
not  damped,  but  a  divergent  oscillation  arises  commonly  known  as  flutter. 
In  this  event,  after  the  initial  disturbance,  each  successive  vibration 
becomes  of  larger  amplitude.  These  oscillations  can  be  dangerously  large 
and  may  lead  to  structural  failure. 

Aerodynamic  forces  on  the  wing  can  arise  from  oscillations  in  vertical 
translation  and  pitching  at  each  wing  span  station,  and  the  coordinates 
for  the  analysis  have  been  selected  for  their  contributions  to  such 
motions.  The  coordinated  degrees  of  freedom  were  considered  under  two 
major  categories:  (1)  symmetric,  where  vertical  and  pitch  motion  of  the 
helicopter,  rigid  flap  of  the  wing  about  its  skewed  hinge,  and  coupled 
flap  bending-torsion  were  considered,  and  where  all  motions  of  the  right 
wing  were  mirrored  symmetrically  by  identical  motion  of  the  left  wing, 

(2)  anti-symmetric  analysis  where  roll  motion  of  the  helicopter,  rigid 
flap  of  the  wing  about  its  skewed  hinge,  and  coupled  flap  bending  torsion 
of  the  wing  were  considered,  and  where  all  motions  of  the  right  wing  were 
accompanied  by  equal  and  opposite  motions  of  the  left  wing.  The  pilot 
adjustable  control  flaps  were  not  considered  to  be  degrees  of  freedom 
because  of  the  rather  rigid  electrical  screw  jack  actuator  between  the 
wing  and  flap.  In  conventional  wings  where  there  are  cable  control 
systems  from  the  flaps  to  the  pilot,  the  flap  pitch  restraint  is  low, 
giving  the  flap  a  low  natural  frequency  and  making  it  an  important  de¬ 
gree  of  freedom. 

The  flutter  analysis  of  the  floating  wing  arrangement  is  further  compli¬ 
cated  by  the  effect  of  sloshing  fuel.  An  exact  analytical  method  for  in¬ 
cluding  fuel  motion  in  the  flutter  equation  is  not  known,  but  several 
approximate  solutions  have  been  published  which  include  the  effective 
masses  and  inertias  for  the  sloshing  fuel.  However,  the  most  directly 
applicable  method  to  the  present  wing  flutter  study  approximates  the 
sloshing  fuel  as  a  change  in  the  chordwise  center  of  gravity. 

This  report  presents  the  flutter  investigations  conducted  for  the  float¬ 
ing  fuel  wings.  Section  II  describes  the  floating  fuel  wing  and  gives 
its  natural  modes  under  various  fuel  conditions.  The  flutter  analysis 
methods  are  presented  in  Section  III,  together  with  detailed  derivations 
in  Appendix  A  and  Appendix  B  for  the  anti-symmetric  and  symmetric  anal¬ 
yses  respectively.  Section  IV  presents  and  discusses  the  flutter  results 
and  views  the  influence  of  fuel  slosh,  and  the  conclusions  of  the  invest¬ 
igation  are  included  in  Section  V.  Numerical  data  including  a  typical 
computer  sample  are  presented  in  Appendix  C. 
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II.  DESCRIPTION  OF  THE  WING 


Each  of  the  wings  is  hinged  freely  about:  an  axis  set  at  a  yaw  angle  of  45° 
with  respect  to  the  longitudinal  axis  of  the  helicopter.  The  hinge,  offset 
angle  is  such  that  a  flap-up  motion  decreases  the  wing  angle  of  attack,  and 
a  flap-down  motion  increases  this  angle,  thus  tending  to  keep  the  wing  at 
a  constant  lift  without  the  helicopter  having  to  carry  the  wing  weight 
Figure  1  illustrates  the  helicopter-wing  system  as  applied  to  an  H-21 . 

The  wing  analyzed  under  this  feasibility  study  has  a  67  ft.  total  wing  span 
and  an  11.2  ft.  chord.  The  flaps  comprise  30%  of  the  chord,  and  are  pilot 
adjustable  by  electrically  actuated  jacks  to  provide  control  of  the  trim 
attitude.  Each  wing  contains  seven  integral  fuel  tanks,  with  sets  of  tanks 
interconnected  as  shown  in  Figure  2.  Each  wing  tank  carries  166.7  gallons 
of  fuel  giving  a  total  capacity  for  both  wings  of  2,343  gallons.  Individual 
wings  weigh  1000  lbs.  empty,  and  8000  lbs.  fully  fueled  The  helicopter  is 
operated  with  the  cabin  empty,  but  with  its  own  fuel  tanks  full  on  take-off 
with  300  gallons,  at  a  gross  weight  of  11,100  lbs.  This  is  below  the  normal 
gross  weight  of  13,500  lbs.  The  helicopter-tank  to  engine-fuel  system 
operates  in  the  normal  manner.  The  helicopter  tank  is  maintained  at  about 
3/4  full  by  a  pumping  system  from  the  wing  tanks  throughout  the  flight.  The 
present  drainage  schedule  is  for  the  tip  tank  pair  to  be  emptied  first,  then 
the  most  inboard  set  of  tanks,  and  finally  the  middle  tank  pair,  so  as  to 
maintain  the  wing  C.G.  near  the  center  of  the  span,  and  avoid  loading  the 
helicopter . 

Based  on  preliminary  structural  and  weight  estimates,  the  elastic  center 
shown  in  Figure  2  is  at  40%  chord,  and  the  fully  fueled  C.G.  is  at  42.4% 
chord.  Detailed  mass,  mass  moment  of  inertia  and  stiffness  distributions 
are  shown  in  Figure  3.  The  pitching  mass  moments  of  inertia  of  the  fuel 
are  taken  to  be  45%  of  their  rigid  value.  The  percentage  figure  is  taken 
from  the  flutter  and  slosh  test  investigation  of  Reference  19,  Figure  56 

Cruise  speed  for  maximum  range  of  the  aircraft  is  80  knots,  the  maximum 
level  flight  speed  is  90  knots  based  on  power  limitation,  and  the  terminal 
dive  speed  is  103  knots  Stall  speeds  are  30  knots  with  wing  empty  and 
70  knots  with  fully  loaded  wing. 

Several  natural  modes  of  the  wing  considered  in  the  flutter  analyses  are 
shown  in  Figure  4.  First  is  the  rigid  body  flapping  mode  of  the  wing  about 
its  hinge.  The  motion  is  restrained  by  an  aerodynamic  spring  resulting  from 
the  change  of  angle  of  attack  with  flap.  For  80  knots  cruising  speed,  the 
natural  frequency  is  1.310  cps  for  the  empty  wing  and  0  342  cps  for  the 
fully  loaded  wing.  Note  that  because  of  the  skewed  hinge,  each  flapwise 
element  of  the  wing  exhibits  both  a  vertical  and  a  pitch  motion  in  this 
rigid  body  flap  mode.  Wing  bending  natural  modes,  the  first  with  flap 
bending  predominant,  the  second  with  torsion  predominant,  are  also  shown 
in  Figure  4.  Other  fuel  loading  conditions  which  vary  the  chordwise  and 
spanwise  C.G.  locations  are  also  given.  The  rigid  body  modes  are  calculated 
directly  in  Appendix  C;  the  flexible  modes  are  calculated  by  an  associated 
matrix  computer  analysis  with  the  numerical  data  given  in  Appendix  C. 
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WING  SCHEMATIC  DIAGRAM 
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III.  METHOD  OF  ANALYSIS 


The  equations  of  motion  for  the  anti  -symmetric  and  symmetric  analyses  are 
determined  from  the  kinetic  and  potential  energies  by  the  use  of  the 
Lagrange  procedure  in  Appendices  A  and  B  Oscillatory  aerodynamic  loads 
are  introduced  through  the  generalized  force  in  each  equation  of  motion 
using  Theodorsen“s  two-dimensional,  incompressible  unsteady  aerodynamics, 
Reference  7,  with  numerical  values  of  these  airloads  taken  from  AF  Report 
4798,  Reference  8  The  equations  of  motion  are.  then  written  in  matrix 
form  to  obtain  a  flutter  determinant.  Symmetric  and  anti-symmetric  deter¬ 
minants  are  presented  in  Figures  5  and  6;  detail  derivations  are  given  in 
Appendices  A  and  B 

It  is  noted  that  the  flutter  derivations  presented  in  the  appendices  con¬ 
sider  the  calculated  wing  modes  as  normal  modes  of  the.  combined  fuselage- 
wing  system.  Consistent  with  this  assumption,  the  products  of  the  gener¬ 
alized  coordinates  between  the  fuselage  and  wing  are  zero  by  the  orthogo¬ 
nality  condition  that  exists  between  normal  modes  However,  wing  mode 
calculations  of  the  present  analysis  were  performed  without  the  fuselage 
properties.  Therefore,  the  calculated  wing  modes  are  not  normal  modes  of 
the  system,  thereby  eliminating  the  assumed  orthogonality  condition  In 
compensating  for  this  inconsistency,  the  flutter  determinants  shown  in 
Figure  5  and  6  include  the  additional  mass  coupling  terms  between  the 
fuselage  and  wing  coordinates. 

In  general,  each  main  diagonal  element  of  a  flutter  determinant  repre¬ 
sents  the  mass  and  inertia  effective  in  a  degree  of  freedom,  the  spring 
restraint  in  the  degree  of  freedom,  and  the  oscillatory  airloads  pro¬ 
duced  by  and  acting  in  the  degree  of  freedom  Each  off-diagonal  element 
represents  the  airload  coupling  existing  between  the  degree  of  freedom 
of  its  determinant  row  position  and  the  degree  of  freedom  of  its  deter¬ 
minant  column  position, 

The  anti-symmetric  flutter  determinant  is  shown  in  Figure  5.  The  degrees 
of  freedom  are 

9  »  fuselage  roll 

«  first  coupled  flap  bending  torsion  mode,  wing 

Hb2  3  second  coupled  flap  bending* torsion  mode,  wing 

The  first  coordinate,  fuselage  roll,  introduces  mass  terms  and  aerodynamic 
loads,  but  does  not  contain  any  spring  terms  because  cf  the  zero  natural 
frequency  of  this  mode  The  two  flexible  modes  are  written  generally,  so 
as  to  accommodate  any  type  of  wing  mode,  whether  it  be  flexible  bending- 
torsion  or  rigid  body  on  the  wing  hinge  Accommodation  of  three  wing 
modes  would  have  been  more  convenient  in  handling  the  present  problem 
(rigid  flap,  first  and  second  coupled  flap  bending  torsion),  but  an  avail¬ 
able  computer  program  with  two  wing  modes  made  the  latter  more  expedient 
It  is  shown  later  that  all  three  pairs  of  these  modes  were  employed,  so  as 
to  cover  the  flutter  coupling  possibilities 


9 


The  symmetric  flutter  determinant  is  shown  in  Figure  6.  The  degrees  of  free¬ 
dom  are: 

Z  -  fuselage  vertical  motion 

'S  =  fuselage  pitch 

Hbf  =  first  coupled  flap  bending-torsion  mode ,  wing 

Hb2  =  second  coupled  flap  bending-torsion  mode,  wing 

The  first  two  coordinates  introduce  mass  and  aerodynamic  terms  but  do  not 
contain  spring  terms  because  of  their  zero  natural  frequencies.  The  two 
flexible  modes  are  used  in  the  same  way  as  in  the  antisymmetric  case. 

Since  it  is  assumed  that  the  oscillations  are  harmonic,  for  example 
Hbi  =  Hbi^i<u5t  and  =  - u e  iA)t,  the  unknown  flutter  frequency  OJ 

appears  in  every  acceleration  term.  Each  equation  is  divided  through  by 
ui2  ,  causing  the  frequency  to  appear  as  a  denominator  of  the  fuselage  ef¬ 
fective  spring  constant  M].ef  f  cJ  Zbl  and  M2eff  b2>  This  results  in  two 
ratios  ‘-°bl/tJJ  and  bl/cA  in  which  the  unknown  flutter  frequency  appears. 

An  additional  complexity  exists,  however,  in  the  numerical  evaluation  of 
the  airloads.  The  airloads  are  all  functions  of  the  Strouhall  number  V/btO 
where  V  =  aircraft  forward  speed,  ;2b  =  airfoil  chord  and  c*>  =  flutter  fre¬ 
quency.  Since  the  flutter  frequency  is  required  in  order  to  determine  the 
airloads,  but  cannot  be  known  until  the  flutter  determinant  is  solved,  it 
is  seen  that  a  trial  and  error  solution  is  necessary. 

A  conventional  method  of  solution  avoids  the  trial  and  error  procedure  by 
artifically  creating  two  unknowns:  7-.  =  (^°bl /<^->)^>  the  squared  ratio  of 
the  first  coupled  bending  frequency  to  the  flutter  frequency,  and 
'S  =  ^  ^bl/&)  d2  *  t^ie  squared  ratio  of  the  two  modal  frequencies.  Since 
each  element  of  the  determinant  is  a  complex  number  by  virtue  of  the  form 
in  which  the  airloads  are  obtained,  the  determinant  can  be  expanded  and 
then  separated  into  one  real  and  one  imaginary  equation.  These  two  equa¬ 
tions  can  then  be  solved  simultaneously  for  the  two  unknowns  7*-  and  . 

Repetitions  of  this  procedure  are  made  for  several  airload  sets  corresponding 
to  several  values  of  the  parameter  V/bt^  .  The  X  and  'f  results,  along  with 
the  corresponding  V/biO  number,  are  converted  to  an  airspeed  V  and  a  ratio 
kJbl/bdb2>  Finally  a  plot  of  V  versus  ration  bl^^  b2  Tnac^e  >  composed 
of  points  for  each  V/b^o  .  The  intersection  of  the  known  ''°bl/^b2  line 
with  the  plotted  curve  determines  the  flutter  speed  V  of  the  aircraft.  The 
slope  of  the  curve  at  the  intersection  is  a  measure  of  the  sensitivity  of 
the  flutter  to  the  magnitude  of  the  modal  natural  frequencies. 
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FIGURE  5 

ANT I -SYMMETRIC  FLUTTER  DETERMINANT 


FIGURE  6 

SYMMETRIC  FLUTTER  DETERMINANT 


IV.  DISCUSSION  OF  RESULTS 


Results  of  the  flutter  analyses  are  summarized  briefly  in  the  table  below,  with  the 
detailed  flutter  curves  shown  later. 
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The  cases  include  full  fuel,  empty,  tip  cells  partially  full  to  produce 
maximum  forward  and  aft  C.G.  in  the  tip  cell,  and  partial  fuel  in  all  the 
cells  so  as  to  p>roduce  the  maximum  forward  and  aft  C.G.  for  the  whole  wing. 
For  each  fuel  configuration  analyses  have  been  performed  for  various 
modpl  combinations,  rigid  and  first  coupled  flap  bending-torsion,  rigid 
and  second  coupled  flap  bending-torsion ,  first  and  second  coupled  flap 
bending-torsion  modes.  The  table  gives  the  critical  flutter  speed  in 
knots  for  the  anti -symmetric  and  symmetric  analyses. 

Cases  1,  2,3  cover  the  fully  fueled  wing,  and  show  a  minimum  flutter 
speed  at  65  knots.  With  rigid  flap  and  first  flexible  flap-torsion  cou¬ 
pled,  flutter  appears  at  67  knots  for  the  anti-symmetric  case  and  90  knots 
for  the  symmetric  case.  Rigid  flap  with  second  bending  torsion  leaves 
the  anti-symmetric  critical  at  65  knots,  but  increases  the  symmetric  to 
160  knots.  Coupling  the  bending-torsion  modes  produces  no  flutter  solution, 
so  this  case  is  stable. 

Cases  4,  5  and  6  review  the  flutter  results  for  0%  fuel  with  each  wing  at 
the  1000  lb  structural  weight.  Again,  as  in  the  previous  cases,  the  modes 
were  paired,  rigid-first  flexible,  rigid-second  flexible  and  first  flexi¬ 
ble  -  second  flexible,  for  the  flutter  investigations.  Both  combinations 
with  rigid  flap  produce  flutter  speeds,  whereas  the  flexible  pair,  as  for 
the  fully  loaded  wing,  indicates  no  flutter  speed.  Flutter  criticals  here 
are  generally  higher  than  for  the  fully  fueled  wing,  reflecting  an  approx¬ 
imately  2  to  1  frequency  increase  in  the  modes. 

Since  the  7000  lb  fuel  load  is  so  large  with  respect  to  the  1000  lb 
structural  weight,  the  problem  of  fuel  slosh  in  the  chordwise  direction 
is  investigated.  While  much  work  has  been  done  in  this  field,  no  com¬ 
pletely  satisfactory  representation  appeared  to  be  available.  As  an 
approximation  of  the  slosh  effect  it  is  assumed  that  the  major  effect  of 
the  fuel  slosh  is  to  move  the  chordwise  center  of  gravity  to  an  extreme 
forward  or  aft  position. 

Cases  7,  8  and  9  assume  that  the  pair  of  cells  nearest  the  wing  tip  are 
drained  until  a  partial  fuel  level  is  reached.  This  portion  of  the 
fuel  is  then  considered  to  be  sloshed  forward  so  as  to  rigidly  fill  the 
forward  portion  of  the  tank,  and  results  in  a  maximum  forward  C.G.  travel. 
The  fuel  quantity  for  this  condition  is  calculated  in  Appendix  C  to  be 
25.8%  of  the  tip  cell,  normal  fuel.  Flutter  analyses,  using  the  same 
modal  combinations  show  stability  with  the  operating  speed  band  for  three 
of  the  six  pairs;  two  cases  show  medium  flutter  speeds,  and  one  case  a 
low  flutter  speed.  The  rigid  first  flexible  symmetric  flutter  speed  cal¬ 
culated  at  55  knots,  is  within  the  operating  speed  and  less  than  the  mini¬ 
mum  for  a  fully  loaded  wing.  Marginal  flutter  conditions  exist  for  other 
modal  combinations,  namely  rigid-second  flexible  symmetrical  and  rigid- 
first  flexible  anti-symmetrical  with  calculated  speeds  near  90  knots. 

Similar  cases,  where  the  fuel  is  sloshed  aft  so  as  to  rigidly  fill  the 
aft  section  of  the  outboard  tanks,  and  move  their  C.G.'s  to  the  most  aft 
position  are  shown  for  identical  modal  pairs  in  cases  10,  11  and  12. 
Appendix  C  calculates  the  fuel  quantity  for  this  condition  to  be  24.2% 
of  the  maximum  tip  cell  load.  In  this  configuration,  a  low  flutter  speed 
exists  in  each  modal  pair  that  includes  the  rigid  flap  mode.  Flutter 
speeds  as  low  as  48  knots  are  calculated  in  this  group  of  cases 
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It  is  clear  from  the  investigation  of  tip  cell  C„G„  variations  that  the 
chordwise  C.G.  does  influence  the  flutter  speed,  lowering  it  from  65 
knots  with  full  fuel  to  a  minimum  of  48  knots,  a  change  of  17  knotSo 
However,  it  is  probable  that  any  configuration  change  which  would  raise 
the  flutter  speeds  of  the  fully  loaded  wing  would  also  be  effective  for 
these  partial  fuel  cases  , 

The  last  six  cases  presume  a  drastic  slosh  condition,  where  in  cases  13, 
14,  and  15,  every  tank  is  partially  loaded  to  25  8%  of  capacity  and  the 
quantity  is  sloshed  fully  forward,  and  to  24.2%  of  capacity  and  this  quan¬ 
tity  is  sloshed  fully  aft,  The  forward  C  G.  produces  flutter  in  only  two 
modal  combinations,  but  the  aft  C  G„  version  produces  low  flutter  speeds 
whenever  the  rigid  mode  is  included  Forward  C  G»  flutter  speeds  appear 
at  95  and  160  knots  when  including  the  second  flexible  mode  with  rigid 
flap.  With  the  adverse  influence  of  the  aft  C„G„ ,  flutter  speeds  occur 
between  35  and  55  knots  for  the  rigid-flexible  mode  pairs. 

Since  the  cruise  speed  of  the  floating  wing  system  is  slated  to  be  80 
knots,  these  flutter  results  are  obviously  unacceptable.  It  will  be 
shown  later  that  these  low  critical  speeds  are  closely  related  to  the 
general  roll  instability  of  the  delta  hinge  arrangement  Some  improve¬ 
ment  could  be  obtained  through  inclusion  of  hydraulic  dampers  at  the 
hinges,  recommended  previously  for  mechanical  instability  control.  The 
latest  design  arrangement  with  a  straight  hinge  parallel  to  the  fuselage 
longitudinal  axis  and  geared  flaps  to  maintain  the  proper  wing  attitude, 
will  probably  not  exhibit  these  undesirable  flutter  characteristics. 
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Flutter  Frequency  Ratio  Plots 


The  results  of  the  flutter  calculations  are  presented  in  more  detail  in 
Figure  7  through  11.  As  described  under  Section  HI,  the  analysis  is  con¬ 
ducted  by  varying  the  airloads  over  a  range  of  the  parameter  V/bw,  and 
solving  for  flutter  speeds  and  modal  frequency  ratio  as  though  the  latter 
were  unknowns.  These  results  are  plotted  as  flutter  speed  versus  frequency 
ratio,  and  when  the  curves  intersect  actual  modal  frequency  ratios,  a  real 
flutter  speed  is  predicted.  The  frequency  ratio  plots  are  presented  in  the 
following  figures: 


Figure 


Fue 1  Condition 
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8 
9 

10 

11 
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1007.  Fuel 
07.  Fuel 

Tip  Cells  Partial  Fuel,  Forward  C.G. 
Tip  Cells  Partial  Fuel,  Aft  C.G. 

All  Cells  Partial  Fuel,  Forward  C.G. 
All  Cells  Partial  Fuel,  Aft  C.G. 


As  an  example,  results  for  the  1007.  fuel  condition.  Cases  1,  2,  and  3  appear 
in  Figure  7.  Three  modal  combinations  are  presented,  l.  Rigid  Flap-First 
Coupled  Mode,  2.  Rigid  Flap-Second  Coupled  Mode,  3.  First  Coupled  Mode-Second 
Coupled  Mode.  The  actual  frequency  ratios  for  these  mode  combinations  are: 


Fuel  Rigid  Flap 

First 

Coupled  Mode.  CPS  Second  Coupled  Mode.  CPS 

100%  0 

4.08  7.48 

Rigid  Flap  = 

0 

»  0 

First  Coupled  Mode 

4.08 

Rigid  Flap 

0 

=  0 

Second  Coupled  Mode 

7.48 

First  Coupled  Mode 
Second  Coupled  Mode 

4.08 

7.48 

=  0.547 

It  is  noted  that  the  rigid  body  flap  mode  of  the  wing  on  its  delta  hinge  is 
taken  to  be  at  zero  frequency  f  !  the  wing  in  a  vacuum.  The  determinant 
solution  contains  airloads  which  provide  the  steady  aerodynamic  spring.  That 
is,  an  up  motion  of  the  wing  on  its  hinge  reduces  the  angle  of  attack,  reduc¬ 
ing  the  steady  lift  and  causing  the  wing  to  drop  on  its  hinge.  Similarly,  a 
down  motion  of  the  wing  increases  the  angle  of  attack  and  the  lift  and  raises 
the  wing  back  towards  its  neutral  position.  If  the  wing  is  analyzed  in  the 
presence  of  this  type  of  air-spring  alone,  the  flapping  natural  frequency  is 
0.16  CPS. 
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In  the  symmetric  solution,  the  zero  frequency  ratio  line  intersects 
the  stability  boundary  curve  in  both  combinations  with  rigid  flap 
showing  a  flutter  speed  of  90  and  160  knots. 

The  symmetric  solution  composed  of  the  first  and  second  coupled  modes 
produces  no  stability  boundary  curve  within  the  normal  range  of  V/bw 
values.  Lower  flutter  speeds  are  obtained  from  the  anti-symmetric 
solution  having  intersections  which  define  flutter  speeds  at  67  and 
65  knots.  Again,  no  stability  boundary  curve  is  produced  by  the 
coupled  mode  combination.  Since  the  V/bw  points  shown  cover  the 
range  of  practical  interest,  additional  points  required  for  defining 
the  stability  boundary  curve  in  some  cases  were  not  calculated. 
Frequency  ratio  curves  for  the  additional  fuel  conditions  are  pre¬ 
sented  using  a  similar  format  in  Figures  8  through  12. 
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Figure  7 


Curves,  100X  Fuel  Condition 
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Figure  8 

Flutter  Curves,  0%  Fuel  Condition 


Airspeed  V  knots 


Figure  9 


Flutter  Curves,  Tip  Cells  Partial  Fuel,  Forward  C.G. 


Figure  10 

Flutter  Curves,  Tip  Cells  Partial  Fuel,  Aft  C.G. 
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Figure  11 


Flutter  Curves,  All  Cells  Partial  Fuel,  Forward  C.G. 


Figure  12 
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Flutter  Curves,  All  Cells  Partial  Fuel,  Aft  C.G. 
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Flutter  Literature  on  Fuel  Slosh 


In  this  investigation,  the  gross  effect  of  fuel  sloshing  is  approxi¬ 
mately  simulated  by  determining  a  most  forward  and  most  aft  extreme 
C.G.  position  for  a  practical  case.  The  maximum  possible  forward  and 
aft  C.G  condition  are  calculated  in  Appendix  C,  and  result  in  the 
following  conditions: 

Partial  Fuel  Cond it  ion s 


3.0"  15.45"  ,  20.8" 


Normal,  100%  Fuel  Most  Fwd  .  C.G.  Most  A.ft  C.G. 

25.8°/  Fuel  24.2 %  Fuel 


While  the  above  method  of  C.G.  shifting  results  in  a  reasonable  solution 
for  the  fuel  slosh  problem  with  the  present  wing,  the  literature  was  re¬ 
viewed  for  other  work  on  the  subject.  Most  investigations  have  been  aimed 
at  determining  slosh  frequencies  and  effective  masses  and  inertias  for  the 
sloshing  fuels,  with  most  of  the  recent  work  directed  at  liquid  fueled 
rocket  applications. 

One  of  the  basic  works  on  the  subject  is  that  of  Lamb,  Reference  11,  who 
showed  that  for  a  rectangular  container  the  slosh  frequency  is  dependent 
on  the  length  of  the  container  and  the  liquid  depth.  If  the  present  wing 
tank  dimensions  are  used  in  the  Lamb  expression,  assuming  only  chordwise 
sloshing  the  slosh  frequency  varies  with  fuel  depth  as  follows: 


0J 


_ l 

2  fr 


1Y _  lanh  (IT  /n) 


where  g,  acceleration  due  to  gravity 
h,  depth  of  fluid 

A  L  ,  length  of  rectangular  container 

_  -  AL 

n  T —  ,  length  to  heignt  ratio 


AL 


CPS 


24 


Frequency  in  CPS 


Direction  of 
Osc illation 


Empty  Full 


Fuel  Quantity  in  Percent  -  % 


An  exact  analytical  method  for  including  fuel  motion  in  the  flutter 
equations  is  not  known,  but  several  approximate  solutions  have  been 
published.  Luskin  and  Lapin,  Reference  12,  attempted  to  include  fuel 
sloshing  in  aircraft  buffeting  analyses  by  a  linearization  in  the 
equations  of  motion.  The  fuel  was  represented  by  a  simple  spring  mass 
system.  Other  works  which  estimate  effective  fuel  properties  are 
References  13  through  19;  Reference  15  also  contains  an  extensive 
bibliography  of  other  works  in  the  field. 

The  most  directly  applicable  reference  to  the  present  wing  flutter 
study  is  Reference  19.  This  is  an  experimental  investigation  of  the 
effect  of  fuel  sloshing  in  the  translation-pitching  flutter  of  a  model 
wing.  Experiments  were  performed  on  a  model  wing  with  internal  fuel 
tanks  restrained  by  springs  to  produce  uncoupled  vertical  and 

pitching  motions.  There  are  three  fuel  cells  arranged  chordwise; 
i.e.,  one  tank  near  the  leading  edge,  the  second  at  the  midsection, 
and  the  third  near  the  trailing  edge.  By  variations  of  the  empty¬ 
ing  sequence,  the  effect  of  chordwise  C.G.  and  percentage  of  fuel 
on  the  flutter  speed  of  the  model  may  be  seen  A  plot  of  the  flutter 
speed  versus  percent  fuel  for  three  emptying  sequences  is  shown  in 
Figure  13. 
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FLUTTER  SPEED  -  FT. /SEC 


FUEL  CELL  EMPTYING  SEQUENCE 
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v.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  possibility  of  flutter  on  the  H-21  Range  Extension  System  wing  with 
the  delta  hinge  arrangement  has  been  investigated  over  a  range  of  fuel 
weight  conditions.  Calculations  have  been  performed  for  symmetric  arid 
anti-symmetric  conditions,  and  for  combinations  of  the  three  lowest  wing 
modes  . 

For  the  wing  with  100%  fuel  load,  flutter  is  predicted  at  65  knots,  within 
the  operating  speed  range.  As  fuel  is  consumed  the  flutter  speed  increases 

to  90  knots  with  0%  fuel.  Compared  to  a  cruise  speed  of  80  knots,  and  a 

maximum  speed  of  90  knots,  both  the  0%  and  100%  fuel  flutter  speeds  are 

unacceptable . 

These  low  critical  flutter  speeds  are  closely  related  to  the  inflight 
static  roll  instability  determined  to  be  characteristic  of  the  delta  hinge 
configuration.  The  critical  flutter  modes  involve  considerable  rigid  body 
motion  of  the  wing  on  its  hinge  at  frequencies  under  0.5  cps.  A  modified 
configuration  currently  under  study  eliminates  the  delta  hinge,  replacing 
it  with  a  straight  hinge  and  a  differential  flap  arrangement  on  each  wing. 

The  latter  provides  an  aerodynamic  spring  similar  to  the  delta  hinge  which 
tends  to  keep  the  wing  in  its  neutral  position,  but  is  designed  so  as  to  be 
statically  and  dynamically  stable.  This  design  change  will  considerably 
modify  the  unusual  flap-pitch  coupling  in  the  oscillatory  airloads,  and 
should  improve  the  flutter  situation. 

A  number  of  partial  fuel  configurations  were  included  in  the  flutter  analy¬ 
sis  of  the. delta  hinge  arrangement.  These  may  be  considered  on  a  relative 
basis,  and  will  point  up  the  fueling  conditions  which  can  significantly  affect 
the  flutter  speeds.  The  influence  of  fuel  slosh  is  investigated  by  first  pre¬ 
suming  partial  fuel  in  the  tip  tanks  with  the  quantity  sloshed  fully  forward  or 
aft  so  as  to  produce  maximum  chordwise  C„G.  travel,  and  then  behaving  as  a  rigid 
mass  in  that  position.  The  forward  C.G.  decreases  the  flutter  speed  to  55  knots, 
the  aft  C.G.  decreases  the  flutter  speed  to  55  knots  relative  to  a  100%  fuel 
flutter  speed  of  65  knots.  A  more  drastic  sloshed  condition  is  also  considered 
in  which  partial  fuel  in  all  tanks  is  sloshed  fully  forward  and  fully  aft.  In 
this  instance  the  forward  C.G.  raises  the  flutter  speed  to  95  knots  and  the  aft 
C.Go  lowers  it  to  35  knots.  Partial  fuel  in  all  the  tanks  is  thus  a  very  un¬ 
desirable  situation.  This  can  be  avoided  by  first  fully  draining  the  tip  tanks, 
then  the  inboard  tanks,  and  finally  the  midspan  tanks.  This  procedure  also  has 
the  advantage  of  tending  to  maintain  the  spanwise  C.G.  near  the  spanwise  center 

of  lift  so  that  loads  on  the  fuselage  of  the  helicopter  are  minimized. 

To  evaluate  the  final  wing  configuration  from  a  flutter  standpoint,  the  safest 

course  to  follow  is  a  build  up  procedure  in  the  initial  flight  testing.  First, 

it  is  recommended  that  a  ground  shake  test  of  the  wings  be  performed  on  the 
ground  with  wings  suspended  in  a  manner  simulating  the  airborne  rigid  body 
spring  rate.  This  will  evaluate  the  wing  frequency  estimates  used  in  the  analy¬ 
sis.  Secondly,  flight  flutter  testing  should  be  conducted  with  sufficient 
strain  gage  and  potentiometer  instrumentation  to  define  the  wing  motion,  and 
a  shaker  should  be  used  to  excite  the  wing.  Initial  flights  would  be  conducted 
with  the  wing  empty,  and  the  measured  data  examined  after  the  flight  for  any 
sign  of  flutter  instability.  Successively  larger  fuel  quantities  would  then 
be  flown,  and  the  data  examined  until  the  full  range  of  fuel  variations  had 
been  explored.  As  an  alternate  to  shaker  excitation,  the  air  jet  method  of 
Reference  20  could  be  used  for  the  inflight  investigations. 
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APPENDIX  A 

ANTI -SYMMETRIC  FLUTTER  ANALYSIS 
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SCHEMATIC  OF  DEGREES  OF  FREEDOM  FOR  ANTI- SYMMETRIC  FLUTTER  ANALYSIS 


Typical  Wing  Section 
Kinematics 
(Looking  Inboard) 
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System  Kinematics 


Wing  vertical  displacement  at  span  station  i, 
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Factoring  out  the  generalized  coordinates, 

*  [  Z  /  Z  Z  2  2  2*,,-  *»v  4„-  ]  Z 

+  [  Z y  Z  C  y  2  *■»*<•  z  Z  ]  Z 

*  [2  ?;  %„■  2  2  X  -4.-  4,.]^  ^  [2  *  24.-  +  2  X  t H 

^  [2  %  %  *  2  4,-  t .  2  2  V  %.  %  +  Z  <.  ^  k  /\ 


The  fuselage  angular  displacement  and  velocity  are: 

<£?  and  <£?  respectively  in  roll. 


Energy  Equations 
Kinetic  Energy 
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where 


M  / Zy  mass  of  each  wing  spanwise  segment 

T'r-  acting  at  station  i. 

mass  per  unit  chord  per  wing  segment . 


Kinetic  Energy  (continued) 

Substituting  the  expression  for 

M »  ‘.•"T  Mut,  r,f> 
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Since  ^  and  are  normal  modes,  then  by  the  principle  of  the 

orthogonality  condition  between  normal  modes,  products  of  the  normal 
coordinates  are  zero.  Thus  the  coefficients  of^/^,j  and  /4^2 

are  zero.  Therefore  the  kinetic  energy  expression  is  reduced 
to 
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Potential  Energy 
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where  » 


^  and  are  effective  spring  rates  of  the  wing 

modes  considered. 
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Lagrange 1 s  Equat  ton 


The  equations  of  motion  are  obtained  by  the  following  operation  due  to 
Lagrange . 


-gL  PlLI 
dt Ibfr  J 


-h 


dv  _  yr_  _  i= 

t 


where: 


=  a  generalized  coordinate 

—  generalized  force  in  the  coordinate 

Application  of  the  above  operation  to  the  kinetic  and  potential  energy 
expressions  yields. 
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And  the  equations  of  motion  are, 


-  <7 

^ 


The  complete  equations  of  motion  require  expressions  for  the  generalized 
forces  A  ,  These  are  derived  in  the  next  section. 
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Generalized  Forces 


The  virtual  work 
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and  the  generalized  force  expression  is 

F-  &JL 

*  Xty 


or: 
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i  (Loct 

Generalized  f-/^  Force 


-»<  Tap 


Generalized  htiz.  Force 
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where : 

£l  —  lift  at  elastic  axis  per  unit  span,  positive  down 

t 

K-  =  moment  at  elastic  axis  per  unit  span,  positive  nose-up 
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Generalized  Forces  as  Oscillatory  Aerodynamic  Loads 
From  AF  Rep  4798  pages  33  and  34, 

L  =  *  ,<* 

M '=  Vi  ]  +  =<  VVV 

where: 

/_  ~  oscillatory  aerodynamic  lift  acting  on  the  wing  per  unit  span. 

M  —  oscillatory  aerodynamic  moment  acting  on  the  wing  about  its  elastic 
axis  per  upit  span. 

<X  -  torsional  displacement  about  the  elastic  axis  of  the  entire  wing, 
positive  for  increasing  angle  of  attack. 

h.  —  vertical  bending  displacement  of  wing,  positive  down. 

(X  —  location  of  the  wing  elastic  axis  measured  from  the  wing  mid-chord, 
as  a  fraction  of  the  semi-chord,  positive  aft  of  mid¬ 
chord  . 

b  —  wing  semi-chord . 


V  —  flutter  frequency. 

P  =  air  mass  density. 

L,  ,  M,  and  A//.  are  complex  aerodynamic  coefficients  tabulated 

;  in  AF  Rept .  No.  4798  for  various  values  of  the  ) 

where  /  =  forward  speed  of  aircraft . 


Rewriting  these  oscillatory  aerodynamic  forces  using  the  notation  of  the 
present  analysis,  the  lift  per  unit  span  becomes. 
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and  the  aerodynamic  moment  about  the  wing  elastic  axis  per  unit  span 
becomes , 
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Complete  Equations  of  Motion 

Substituting  the  oscillatory  aerodynamic  coefficients  into  the  generalized 
force  expressions^ then  tfye  equations  of  motion  become, 
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Collecting  the  coefficients  of  the  generalized  coordinates  and  assuming  a 
harmonic  solution  of  the  form: 


t  =  / 


'  ejJ  t 


such  that: 


eCurt  --cS-Q- 


/ 


6 


then  upon  substitution, the  flutter  equations  may  be  written  in  Matrix  form 
as  shown  in  Figure  4 


where: 


A  =  M L 

4  = 


2.GFF 


- 


natural  frequency  of  the  first  wing  normal  mode 
in  rad/sec. 


<■0, 


?. 


natural  frequency  of  the  second  wing  normal  mode 
in  rad/sec. 
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Solution  of  the  Antisymmetric  Flutter  Matrix 

For  a  solution  to  exist  the  determinant  of  the  flutter  matrix  coefficients 
must  be  zero.  Since  the  coefficients  are  complex,  expansion  of  the  determinant 
results  in  two  equations  with  reai  coefficients,  and  only  one  unknown,  the 
flutter  frequency  QQ  .  Thi,s  means  that  one  equation  is  extraneous .  To  circum¬ 
vent  this,  assume  another  variable,  which, in  this  analysis  is  the  ratio  of  the 
frequencies  of  the  two  wing  modes.  The  two  equations  may  then  be  solved  simul¬ 
taneously  for  tw'o  unknowns.  Repeated  solutions  of  t;he  flutter  determinant  in 
this  manner  must  be  performed  for  various  values  of  the  reduced  frequency  V/b(A . 
The  solutions  of  the  equations  may  then  be  converted  to  a  plot  of  the  forward 
velocity,  (/  against  the  corresponding  frequency  ratio,  from  which  the  flutter 
speed  may  be  obtained  by  the  intersection  of  the  actual  frequency  ratio  with 
the  curve. 

The  steps  for  the  solution  of  the  flutter  determinant  are  as  follows.  Let 
Roman  Numerals  represent  the  coefficients  of  the  determinant  (remembering 
that  these  are  complex) . 
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^  -  mJ&* 


IV 
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Now,  divide  the  second  row  by  MleFf.  and  the  third  row  by  Fjtp  >  and 

represent  the  resulting  terms  by  A,  3jC,  . .  which  are  still  complex. 

Then  making  the  following  substitutions. 


and 
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_  (“>h, 


the  determinant  becomes 
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£  -  A 


H  t 
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F 


A 


=  O 


Solution  of  this  complex  determinant  for  A  and  S  is  programmed  (Program 
No.  1120)  on  the  IBM  650  computer,  and  the  output  is  as  follows. 
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Antisymmetric  Flutter  Program  -  Output  Format  (Program  No.  1120) 


Ar 


Br 


Cr 
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Dr 
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Er 

El 
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GI 


Hr 

hi 


Kr 

Ki 


vknot  s 


Zeta 


Lambda 


where ,  A ,  B 


K  =  coefficients  of  the  reduced  determinant. 

Subscript  R  denotes  the  real  part; 
subscript  I  denotes  the  imaginary  part . 

M,  N,  P  =  coefficients  (real)  of  the  equation  in  /V.,e.i. 

~K  M  +  A  A/  +-  P  O 


^knots 


Zeta 


=  forward  speed  in  knots; 


1/ 


hu>  /aT 


x  Conversion  factor  to 
knot  s . 


Lambda 
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APPENDIX  B 

SYMMETRIC  FLUTTER  ANALYSIS 
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Typical  Wing  Section 
Kinematics 

(Looking  Inboard,  Left  Wing) 


Symmetric  Mode 


System  Kinematics 


The  wing  vertical  displacement  at  Station  i. 


2  ~  s  y  + 


The  fuselage  C.G.  displacements: 
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and  the  velocities  are: 
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Since,  2,  ^  and  A^,  are  normal  modes,  then  by  the  principle  of 

orthogonality,  coefficients  of  the  products  of  the  generalized  coordinates 
must  be  zero.  Applying  this  principal  on  the  last  equation  yields. 
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Symmetric  Modes 


Energy  Equations 


Kinetic  energy 
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Mw.  -  J  <y,  mass  of  each  wing  spanwise  segment 

^  acting  at  Station  i 

mass  per  unit  chord  per  wing  segment 


Substituting  the  expressions  for 
becomes , 


the  kinetic  energy 
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Collecting  coefficients  of  the  generalized  coordinates  gives, 
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Symmetric  Modes 


Kinetic  Energy  (continued) 

The  kinetic  energy  may  be  rewritten  thus, 

T  -  z\f  B  *  C  +  D~fx\ 


Potential  Energy 


x  "V^ 


where: 


effective  spring  constant  of  each  wing  mode. 


Lagrange's  Equation 

The  mechanical  equations  of  motion  are  obtained  by  the  following  operation 
(due  to  Lagrange)  on  the  kinetic  and  potential  energy  expressions. 


where : 


d 

dt 


dT 


+ 


bv  _  az 


generalized  coordinate 


generalized  force  in  the  coordinate 


f 


These  result  in  mechanical  equations  of  motion  with  the  same  number  as  there 
are  generalized  coordinates. 


Substitution  of  the  kinetic  and  potential  energy  expressions  into  the  dif¬ 
ferential  operation  above  yields  the  following  equations  of  motion. 
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Symmetric  Modes 


Lagrange's  Equation  (continued) 
Equations  of  action, 


A. 


.  OO 


Generalized  Forces 


To  complete  the  equations  of  motion,  the  generalized  forces  must  be  found 
These  may  be  obtained  from  the  virtual  work.  Thus, 


£  = 

aw 

If  = 

virtual  work 

rr 

•  jp.  and 

are  as  defined  above. 

increment  of 

where: 


The  virtual  work  is, 

0  ,  TlP 

I 


'  C 

Lc  4 ^  +  Sb-  S  $  X  v 


and  the  generalized  forces  are  then 

w,  tv/5 


c=/j 

n  t/p 


i=/,  et&T 

n,  t/p 


n  7y/° 


t'-i 


Generalized 

a 

force 

Generalized 

y 

force 

Generalized 

fit, 

force 

Generalized 

Hbz 

force 

46 


Symmetric  Modes 

Generalized  Forces  (continued) 
where : 

^■i  -  lift  at  elastic  axis  per  unit  span,  positive  down 

/ 

moment  about  elastic  axis  per  unit  span,  positive  nose-up. 

Using  AF  Report  No.  4798,  the  generalized  forces  may  be  rewritten  as  oscil¬ 
latory  aerodynamic  loads . 


Generalized  Forces  as  Oscillatory  Aerodynamic  Loads 
From  Reference  (8)  pages  33  and  34, 

L  -  7T  pi,  of  {~t~  \Li,  ]  +  (±+a)  ]  } 


M  -  La  (-£+a)]  ]} 


where: 


Z  -  oscillatory  aerodynamic  lift  acting  on  the  wing  per  unit  span 

M-  oscillatory  aerodynamic  moment  acting  on  the  wing  about  its 

elastic  axis  per  unit  span 

o(  —  torsional  displacement  about  the  elastic  axis  of  the  entire 

wing,  positive  for  increasing  angle  of  attack 


In 


vertical  bending  displacement  of  wing,  positive  down 


CL  —  location  of  the  wing  elastic  axis  measured  from  the  wing  mid¬ 
chord,  as  a  fraction  of  the  semi-chord, 
positive  aft  of  mid-chord 

b  -  wing  semi-chord 


~  flutter  frequency 
P  —  air  mass  density 


Z  >  > 


K 


and  are  complex  aerodynamic  coefficients  tabulated 

^  in  AF  Rept .  No.  4798  for  various  values  of  the 
where  =  forward  speed  of  aircraft. 
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Symmetric  Modes 


Generalized  Forces  as  Oscillatory  Aerodynamic  Loads  (continued) 

Rewriting  these  oscillatory  aerodynamic  forces  using  the  notation  of  the 
present  analysis,  the  lift  per  unit  span  becomes, 

L  -  ff/°jb3coZ [  ± z**ttf**.  y-a)J  j 

and  the  aerodynamic  moment  about  the  wing  elastic  axis  per  unit  span  is, 


V  ‘  2i 


Complete  Equations  of  Motion 

Substituting  the  expressions  for  the  generalized  forces  in  the  equations  of 
motion  and  writing  out  A,  B,  C  and  D  in  full,  the  complete  equations  of  motion 
are  obtained.  Thus, 


Equation 


r-  1.T/P  "iTir 

J  1  ^  =  y  rrp  ban  *  -  gx ±  ;/7ih 1 d**  p.,1 

/  */  JZooT'  — 1  .  ^7  ■  y  -  ^ 

4  -*  i  R.oot  4 


2f  Equation 


.  ,  n,  T  i  P 

ri-i  *  z^,s*i?= 

L  i«lj  J 


^  p b  6oj  ^  -sy+  //^  p/ £.x 

1  raf] 

S  /  b3 td"1' j  SY+^i/? J.  j"^  j  1 

V  (r  ^  ^  ^ C  i + *)]  j y  ‘ 
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Symmetric  Modes 

Complete  Equations  of  Motion  (continued) 
Equat  ion 


y)j  tip 

i^/jKDO  7-  '*  ' 

»?,  TVg 

7>^V4, (  |i'4-  /.4<ri^] 


■/JooT 


rdyt. 


^  Equation 


TIP 


3  ^  t,  -  -t,  C  *■ 2  %  .  ]  4  *-  4 a= 

1  *i)  eooT  c  A  2^  ; 


e»o 

*,  HP 


^Firptfur^ _  j'  g-sy  *  g*^//**  -  t-A(-±  +  aj] 


(car 


+ 


<*+4  HOty*  4,  i  +«) +Lci «)1 


Collecting  the  coefficients  of  the  generalized  coordinates  and  assuming 
harmonic  solution  of  the  form 


/' 


so  that 
=  “"> 


then  upon  substitution,  the  equations  of  motion  or  flutter  equations  may 
be  written  in  matrix  form  as  in  Figure  5; 


where: 


■^/V>  *  i, 

F 

W-  A. 


cj 


Z£Ff 

GJ  ~  natural  frequency  of  the  first  wing  normal  mode 

b, 

CJ  natural  frequency  of  the  second  wing  normal  mode 
b 
i 

(JJW’  __  ___ 

V  l-ljftooT 

'i  =  n^,Ar  *A* +  ^ 

i— I,  Root 
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Symmetric  Modes 


Solution  of  the  Symmetric  Flutter 

The  solution  of  the  Symmetric  Flutter  determinant  is  identical  to  the 
solution  of  the  Antisymmetric  Case  with  one  exception:  instead  of  a  3  x  3 
determinant,  the  Symmetric  Flutter  determinant  is  a  4  x  4.  This  must  be 
reduced  to  a  3  x  3  and  may  then  be  solved  in  a  manner  identical  to  the 
Antisymmetric  Solution.  This  reduction  of  the  4x4toa3x3  deter¬ 
minant  is  done  as  follows. 

Represent  the  elements  of  the  determinant  by  Roman  Numerals. 


1 1  nr 

[.  E  ; 

I 

* 

xi 

m 

ML 

IX 

x  j 

XL 

M  ~ 

HE 

w 

XX  1 

szi:  -  n  tej 

Divide  each  element  of  each  column  by  the  first  (top)  element  of  the  column 
This  reduces  the  first  row  elements  to  unity. 


/ 

/ 

/ 

/ 

YL 

XL 

m 

x 

in 

X 

1 

X 

X 

XX 

XEL  I  Utf 

m 

XL 

XX 

m  ez:IW 

X 

m 

XX 

w 

XVLL  C., 

TL 

XX 

u 

XX  X1 
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Symmetric  Modes 


Solution  of  the  Symmetric  Flutter  Determinant  (continued) 


Subtracting  the  first  column  from  the  other  three  columns  reduces  the  first 
row  elements  to  zero,  except  the  first,  which  is  still  unity.  Therefore, 
by  Cramer's  rule  the  4x4  determinant  reduces  to  the  minor  of  the  first- 
row-column  element,  1. 


m 

m 

SID  ST 

IX 

m 

n 

EZ  nr 

ez: 

rn 

XI 

X 

XU  X  K*  l^)z 

XU 

X 

nr 

nr 

ZEL  XL  EZ :  W 

3 L 

rr 

xy 

m 

m  M 

m  m 

m 

rr 

m.  nr 

eh 

hi 

311“/ 

The  above  determinant  has  a  form  identical  to  the  Antisymmetric  case;  hence 
its  solution  also.To  put  the  reduced  matrix  into  its  final  form,  divide  the 
second  column  by  the  coefficient  of  /cO  )3~  and  the  third  column  by 

the  coefficient  of  (  CO ^ ^  CO>y~  •  Also  make  the  following  substitution: 


'2. 

A  and  f  so  that 

WJ 

After  further  multiplying  the  third  row  by  ^  ,  the 

becomes 


Kf-  2A 
W  /  ~  S 

determinant  finally 
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Symmetric  Modes 


Solution  of  the  Symmetric  Flutter  Determinant  (continued) 


M  XL 


nr 

XT 

XL 

X 

rx 

XL 

“XX 

W 

hi  xr 

r- 


H  . 


M, 

*rf 


M 

V 


2m- 


XEL- 


X>1X 


XT 


XjlE  i  TV 
IX  _j 


/ 

'  Mz 

ere 

H. 


LX- 


w 


X2 


:xl  x 

XL 


x  x 

XT 


m  x 

XT 


f-a 


This  may  also  be  written  in  an  abbreviated  form  by  representing  the  elements 
by  capital  letters.  It  should  be  noted  that  the  elements  of  the  determinant 
are  complex. 


A 


D 


3 


I  -  A 


//y 


-  A. 


=  O 


The  solution  of  the  Symmetric  Flutter  determinant  consists  of  (1)  reduction 
of  a  4  x  4  determinant  to  a  3  x  3  and  (2)  simultaneous  solution  of  the  two 
equations  resulting  from  the  expansion  of  the  reduced  determiant  for  the 
pseudo  unknown  frequency  ratio  ~f  and  7^-  ,  is  programmed  on  the  IBM  650 

computer.  The  program  is  designated  No.  1178.  The  output  format  is  the 
same  as  in  the  Antisymmetric  case,  which  was  shown  previously. 
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APPENDIX  C 


NUMERICAL  DATA 


TABLE  I 


H-21  RANGE  EXTENSION  HELICOPTER,  BASIC  NUMERICAL  DATA 


Symbols 

Dimensions 

Value 

Description 

P 

lb-sec^/ in^ 

0.1147  x  10"6 

Air  mass  density  @  sea  level 

a' 

in 

-13.4 

Wing  midchord  to  E.A.  positive 
aft  of  midchord 

a 

- 

.199 

a'/b 

b 

in 

67.4 

Wing  semichord 

V 

ft/sec 

135.0 

Forward  speed  of  aircraft 

S 

in 

61 .0 

Aircraft  C.G.  aft  of  wing  E.A. 

M/2 

lb-sec^/ in 

14.4 

Half  of  fuselage  mass 

XW2 

lb-in-sec^ 

30,000 

Half  of  fuselage  moment  of 
inertia  in  roll 

Xfp/2 

lb-in-sec^ 

404,000 

Half  of  fuselage  moment  of 
inertia  in  pitch 

X 

deg. 

45.0 

Wing  hinge  offset  angle 

in-lb/rad 

3.09  x  106 

Wing  aero,  spring  rate  at 

80  knots  (135  FPS) 
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DETERMINATION  OF  FLEXIBLE  AND  RIGID  WING  MODES 


The  frequencies  and  corresponding  wing  flexible  mode  shapes  are  calculated  using 
the  method  of  associated  matrices  programmed  in  the  IBM  650  Computer  For 
the  mode  shape  calculations*  the  wing  is  represented  as  a  series  of  discrete 
lumped  masses  interconnected  by  massless  elastic  members  which  have  both 
bending  and  torsional  rigidity.  Boundary  conditions  are  for  those  of  a  pinned- 
free  beam  and  the  resulting  modes  are  coupled  bending  and  torsion  which  are 
also  normal  wing  modes.  The  first  flexible  mode  is  predominantly  flap  bending 
and  the  second  mode  is  usually  torsion  predominant.  Inputs  are  prepared  in  the 
sample  input  sheet  labeled  Table  (  II  ) . 

Due  to  the  wing  hinge,  a  rigid  mode  of  the  wing  exists  during  forward  flight, 
and  which  consists  in  flapwise  rotation  about  the  hinge.  Determination  of  the 
wing  rigid  mode  is  made  with  the  aid  of  projective  geometry. 

Referring  to  the  figure  in  the  next  page  which  shows  the  deflected  wing  as 
solid  lines  and  the  equilibrium  position  by  dotted  lines,  the  following  are 
obtained: 

For  small  wing  flapwise  deflections, 
d3c3-  33c3=  d^ 

but 


dlCl  =  ai°itan  ^  o 
Also: 

030.3'  =  c2c2'  =  c2b2sin  (V  w  =  c2b2sinO<w  =  Cjt^sin  Y  w  cos 
and 

d3d3'  =  d2d2  =  d2a2 '  sin°^w  =  dlal  sin°^w  cos  Y0 
Now 

Sin  0...  =  c,c,'  -  d_d_  or  after  substitution 

'  bli  3  3 _ 3  3 

d3c3 

Sin  ^bli  =  Uibi  -  d!a!)  Sin*w  COS^  o 
a^o^  tan 

=  YiSinO(w  Sin  Y  o 

al°l 

Sin^bli=  Sin  ^  W  Sin  Y  o 
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RIGID  WING  FLAP  GEOMETRY 


Primed  letters  denote  projections  of  final 
position  in  the  equilibrium  plane 


and  since  angles  are  small 


3  =  a  sin  v 

'bli  w  6  o 

The  vertical  deflections  are 

=  y  cosy 

i  0  3 


bli 


COS  ^ c  bL  Tip  R  bl  Tip 


Where  R  is  the  length  of  the  elastic  axis  from  wing  tip  to  the  hinge 
unit  tip  deflection. 


7-  bli 


=  n 


Also,  since 

?bi  '  1  '  <(“),)»(  . 

tl  en: 

OC  =  1 

w  R  COS'tf 

to  that 

=  Tan  X  q>  constant  along  the  span. 
R 

Now  for: 

R  =  303" 

2. 


bli 
^  bli 


yo  =  45° 

=  y i  ( .0033) 

=  .0033  Rad.,  constant  along  the  span 
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For  a 


TABLE  II 
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TABLE  II  (Continued) 

BASIC  NUMERICAL  DATA  (Continuation) 


SECTION  ELASTIC  PROPERTIES  OF  THE  WING 


Station 

—  in, 

in 

EI-lb-in2 

GJ~lb-rin2 

Tip  1 

1 

303 

2. 50x1 09 

5. 38x1 09 

2 

mm 

258 

3 

28 

213 

4 

28 

185 

5 

45 

157 

6 

45 

112 

7 

45 

67 

Root  8 

22 

22 

y 

Schematic  Diagram  of  the  Wing  for  Mode-Shape  Calculation 


The  derivation  and  inputs  for  the  mode-shape  calculation  were  previously 
obtained  using  symbols  different  from  the  Flutter  nomenclature.  In  the 
sketch  below  the  mode-shape  calculation  symbols  are  shown,  including  the 
sign  convention. 
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When  the  aircraft  is  airborne,  the  wing  can  execute  rigid  flapping  about 
the  hinge  at  a  certain  frequency.  The  spring  rate  is  due  to  aerodynamLc 
loads  on  the  wing.  To  determine  this  spring  rate,  consider  the  following. 


Wing  Hinge  / 


Then  from  the  "strip"  theory,  the  steady  flapping  moment  about  the  wing 
hinge  axis  is , 

m  =7?/°<z.  at)i'*0ycos%4 

*  o 

from  the  calculation  of  the  rigid  mode  shape, 

^  S/ai'Y  -  substituting  this  above,  and  simplifying, 

w  °  e 

M*w  =  ^  S/fl  yo  Cos'd' D yS  c/y  -  upon  integration 

=  4  IZ70 

This  may  also  be  rewritten  as: 

M«W  ‘  Kfjx*  where 

“  the  effective  aerodynamic  spring.  (This  has  been  calculated  in 
Reference  (6)  at  sea  level  and  80  knots  to  be: 

=  3.09  x  10^  in-lb/rad.) 
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DETERMINATION  OF  MAXIMUM  CHORDWISE  DISPLACEMENT  OF 
WING  C  G .  DUE  TO  CHANGE  IN  FUEL  LEVEL 


Aft  -  Displacement 

The  H-21  Range  Extension  Wing  has  six  fuel  cells  arranged  along  each  semi¬ 
span.  The  procedure  for  transferring  the  wing  fuel  to  the  main  tank  is  as 
follows: 

(1)  The  first  two  fuel  cells  nearest  the  wing  tip  are  simultaneously 
transferred  to  the  main  tank. 

(2)  When  the  first  two  cells  are  empty,  the  most  inboard  two  cells 
are  emptied.  Then  the  midspan  tanks  are  emptied  last. 

To  determine  the  percentage  fuel  remaining  in  a  cell  that  yields  maximum 
displacement  of  the  wing  section  C.G.  in  the  chordwise  direction,  assume 
that  any  fuel  in  the  cell  which  has  a  normally  horizontal  free  surface 
occupies  the  rear  end  of  the  cell.  In  this  condition  the  fuel  would  have 
a  vertical  free  surface. 

Considering  the  wing  to  be  horizontal  and  the  cell  cross  section  to  be  rec¬ 
tangular,  the  following  sketch  shows  the  geometry  of  the  wing-fuel  cell 


Where:  Ms  =  Mass  of  Unit  Length  of  Wing  Structure 

Mf  =  Mass  of  Fuel  per  Unit  Length  (spanwise) 
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Aft 


From  the  Sketch,  the  Fuel  Mass  is 

M'f  «  2Z  &  =  t( L-  X)  S 

Where,  &  =  Fuel  Mass  per  unit  area 
From  which , 

X  -  L(l-^) 

Also,  it  is  easily  shown  that 
Xf  =  L(l-^) 

The  C.G.  of  the  entire  section  then  becomes, 

X  =  MfXf  a£ter  substituting  Mf  and  Xf 

.  2 


Ms+  and  simplifying, 


X  = 


_  MSXS  + 


Ms  +  2  LS 

will  be  a  maximum  or  a  minimum  when  d  X-/ d  —O  ; 

_  (Ms  +  2L  S)ft  2S-  *l/t}  -//%  ^Vzt)  A.S 


(Ats  ^  zxS)* 


Considering  only  the  numerator  and  solving  for 
obtained . 


the  following  is 


*=-f 


Afs 


~/± 


/ 


/ 


2  t  U-*s) 

Ms  /$ 


] 


And  the  percentage  of  fuel  remaining  is 

-i.  x  ioo 


It  must  be  noted  that  only  positive  values  of  2  is  valid. 


Positive  ^  is  obviously  a  maximum. 

Substituting  this  into  the  expression  for  pi  gives  the  wing  section  maximum 
chordwise  C.G.  travel  Aft. 

The  distance  of  the  wing  C.G.  from  the  E.A.  is  therefore: 


X.  -  X  -  34.2 

S-AtC 
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Forward  Displacement 


To  determine  the  maximum  chordwise  wing  section  C.G.  displacement,  we  pro¬ 
ceed  in  a  similar  manner.  This  time,  assume  the  fuel  to  occupy  forward  end 
of  the  cell  thus: 


Similarly,  then 

M  f  -  2  L  £  =  t  XS 

Which  yields 


And 


The  C.G.  of  the  entire  section  is  then, 

-t-  2.  L  S 


Differentiating  with  respect  to  -g  for  maximum-minimum  and  equating  to 
zero,  we  have 


dx_ 

c/i 


( +  ?  z  Y2t)LS  _  Q 

(M5  f  Zr  /.&)*- 
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Forward 


Considering  only  the  numerator  and  solving  for  2 
obtained: 


ZXsiZ 

Ms/$ 


Again  only  a  positive  value  is  valid. 
This  positive  2  is  a  minimum. 


the  following  is 


Substituting  into  the  expression  for  X  yields  the  minimum  C.G.  chordwise 
displacement,  and  this  distance  with  respect  to  the  elastic  axis  is: 


X  =  X  _  3  4.2 

c-&  -  e.a 


Calculation  of  Maximum  C.G.  Displacements 

Substituting  the  following  parameters  into  the  previous  equations  yield 
the  percent  fuel  and  its  C.G.  location  for  maximum  aft  and  forward  dis¬ 
placement  . 


t  -  18.42" 

L  =  72.6" 

Most  Aft  C.G. 
lx  100  =  24.2% 


Ms  =  .844  lb-sec^/in 
Xs  =  38.2  in 


£  =  5.04  =*A2 


X  =-15.45"  (or  15.45" 
forward  of 
E.A.) 


*  =  20.8"  aft  of  E.A. 

C.  &.A. 


Most  Forward  C.G. 
x  100  =  25.8% 
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TABLE  IV 

ANTISYMMETRIC  FLUTTER  ANALYSIS 
Rigid  and  First  Flap  Bending  Mode 


Mechanical  Inputs 

100%  Full  -  H-21  Range  Extension  Wing  #1 


Program  No.  1120 


Case  4 


TABLE  V 

SYMMETRIC  FLUTTER  ANALYSIS 

Rigid  and  First  Flap  Bending  Mode 
Mechanical  Inpu t s 

100%  Fuel  -  H-21  Range  Extension  Wing  #1 


Program  No.  1178 
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TABLE  VI 


JOB  #4485 


V/BW 

8000000  51 


V/BW 

9000000  51 


TYPICAL  OUTPUT  SHEET  -  FLUTTER  ANALYSIS 
ANTISYMMETRIC  FLUTTER  ANALYSIS 
RIGID  AND  FIRST  FLAP  BENDING  MODES 


CASE  IAS 

Oct.  12, 

1960 

A 

B 

C 

73535170 

56 

-19454193 

53 

85616989 

53 

-23706330 

56 

-76618636 

52 

-10208978 

53 

D 

E 

F 

-90091250 

52 

49577181 

50 

22235164 

51 

-65250324 

53 

-22867608 

50 

-23009971 

50 

G 

H 

K 

-13557957 

53 

-43468827 

50 

27799010 

51 

-68238753 

53 

-17420489 

50 

-71499957 

50 

M 

N 

P 

-36663790 

61 

59048583 

62 

-47293843 

62 

V 

ZETA 

LAMDA 

19852195 

52 

55677156 

51 

15260126 

52 

84349603 

52 

20369645 

51 

84529655 

50 

A 

B 

C 

73011272 

56 

-24751878 

53 

10978946 

54 

-27335358 

56 

-84664698 

52 

-13523939 

53 

D 

E 

F 

-10451124 

53 

35802261 

50 

28520329 

51 

-75239016 

53 

-25378434 

50 

-31256640 

50 

G 

H 

K. 

-15065998 

53 

-54792397 

50 

32852389 

51 

-78684923 

53 

-19482675 

50 

-85332875 

50 

M 

N 

P 

-17151336 

62 

10816808 

63 

-76407132 

62 

V 

ZETA 

LAMDA 

37214423 

52 

18072757 

51 

54961385 

51 

96906073 

52 

34569478 

51 

81054705 
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